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We study the effect of enforcing exact conservation of charges in statistical models of particle 
prodnction for systems as large as those relevant to relativistic heavy ion collisions. By using a 
numerical method developed for small systems, we have been able to approach the large volume limit 
keeping the exact canonical treatment of all relevant charges, namely baryon number, strangeness 
and electric charge. Hence, we hereby give the information needed in a hadron gas model whether 
the canonical treatment is necessary or not in actual cases. Comparison between calculations and 
experimental particle multiplicities is shown. Also, a discussion on relative strangeness chemical 
equilibrium is given. 
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I. INTRODUCTION 

Over the era of relativistic nuclear collisions, the 
statistical-thermal models have been widely used in the 
analysis of particle production. During the past few 
years, these models have been successful in describing 
particle multiplicities in high energy nuclear reactions at 
the stage where inelastic collisions between hadrons cease 
(chemical freeze-out), see e.g. Qiiii 

The easiest to handle, thence the most used, statisti¬ 
cal model calculational framework is based on the grand 
canonical (GC) ensemble. In the GC approach entropy 
is maximized using the constraints of conserved ensem¬ 
ble averages of energy and charges. This allows the net 
charges to fluctuate from sample to sample even though 
the actual charges brought into the physical reactions 
were exactly same every time. In the large volume and 
energy limit, the error in such an assumption due to 
charge fluctuations is negligible. 

On the other hand, when the involved volumes are 
small, like in e“'"e“, pp, pA reactions or in peripheral AA 
reactions, the charge fluctuations allowed in GG formal¬ 
ism give a large theoretical error. So, the GG ensemble 
turns out to be inadequate for the analysis of experi¬ 
mental results and the statistical ensemble to be used is 
rather the canonical one, enforcing exact conservation of 
relevant charges. In this case entropy is maximized using 
the constraint of conserved ensemble average of energy, 
but charges are not allowed to fluctuate. 

Statistical models based on the canonical ensemble 
have been able to reproduce particle multiplicities even 
in elementary collisions, although in those there is clearly 
no room for kinetic thermalization i u Also, the 
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strangeness enhancement with increasing system size in 
p-A reactions can be explained by a canonical effect 0. 

In this work, we investigate whether and to what ex¬ 
tent the canonical treatment is relevant for actual heavy 
ion collisions. The main difficulty of canonical calcula¬ 
tions is related to the involved large values of baryon 
number and electric charge, which make the numerical 
computation of canonical partition functions quite hard. 
We present here an efficient numerical method to carry 
out such calculations which enabled us to implement ex¬ 
act charge conservation up to baryon number 0(100) 
whereas with old methods it was only possible to reach 
B ^ 20 0 with very long computation times. 

After deriving the needed expressions of mean particle 
numbers, we present a systematical study of the chemi¬ 
cal factors appearing in the canonical statistical model. 
Along with some comparison between grand-canonical 
and canonical calculations, we show a comparison with 
experimental results for peripheral to central A-A reac¬ 
tions at AGS and SPS energies. 

In section nvi we briefly address the role of the so- 
called relative strangeness chemical equilibrium, usually 
parametrized with 75, in canonical statistical-thermal 
models. 


II. THE ANALYTICAL DEVELOPMENT 


In nuclear physics, the problem of calculating the rel¬ 
ativistic canonical partition function has been handled 
by several different methods, see e.g. i HI HI HI 
Here, we employ a very general group theoretical method 
- applicable to any internal symmetry represented by 
a semi-simple Lie algebra - first introduced by Gerulus 

El HI HI 


By denoting the set of conserved quantum numbers by 
{Ci}, the canonical partition function can be ob¬ 

tained from the usual grand-canonical partition function 
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Zgc by using a projecting operator onto the conserved 
quantum numbers. In the case of N internal symmetries 
of type the projection takes the form: 


Z{c.}{T,V) 



Zgc{T,V,{\c,}). 


( 1 ) 

where (/>i e [0, 27r) is a 17(1) group parameter and a Wick- 
rotated fugacity factor Ac^ = is introduced, for every 
charge Ci. 

In order to calculate Zgg in tbe above integral, we will 
assume the Boltzmann statistics for all hadron species. 
The error involved in such an approximation for final par¬ 


ticle multiplicities, when resonance decays is taken into 
account, is some percent in case of pions and much less for 
all other hadrons. Quantum statistics may be used, but 
requires somewhat heavier calculations and disregarding 
it does not essentially affect any forthcoming argument. 

In heavy ion reactions, the relevant set of con¬ 
served charges is {Q} = B,S,Q, namely baryon 
number, strangeness and electric charge. Denoting 
the partition function of hadrons carrying none of 
the previous charges by Zq, and writing Zgg by us¬ 
ing the one-particle partition function zj = [{2Ji + 
l)F/(27r)3]/d3pexp[-v^ -|-m^/T] for each hadron i, 
we find 


r2Tr 


Zb,s,q{T,V) = Zo 


(27r)3 


d^B e 


p‘ 2 'k /‘ 27 r 

^-iS4>s I 


dfjjQ e 


-iQ4>Q 


X exp < g^{Bi(j>B+Si^s+Qi4>Q) ^-KBi4>B+Si4>s+Qi^Q) 


( 2 ) 


The direct numerical computation of the triple integral 
above is a formidable task for B >5. The integral can be 
turned into a sum over many indices of modified Bessel 
functions In{x), yet its evaluation is very time consum¬ 
ing and becomes impractical with baryon number larger 
than ^ 10 0. Therefore, we have chosen another ap¬ 
proach, which has been applied in elementary collisions 
[3, Il7| . From m one finds the integrand to be violently 
oscillating for large net quantum numbers, so it is worth 


trying to eliminate analytically the source of strongest 
oscillation; in heavy ion reactions, the baryon number is 
always the largest of the set B, S,Q, so it is certainly the 
most beneficial to eliminate the integration over B. In 
fact, this can be done by taking advantage of a special 
feature of baryon number, that is no elementary hadron 
exists with \Bj \ > I. We first rewrite equation © in the 
form 


-I p27r p27T 

Zb,s,q{T.V) = Zo-^ / d()>Be-‘^'^^ / / d^g 

(27r) 7o JO Jo 


X exp 


exp 



_b,'b 


_M.M 


where the two summations run over baryons and mesons, 
respectively. As \Bi\ = 1 for all of the baryons, we can 
write the baryon summation as: 


^1 gi{Bi4'B+Si4>s+Qi'pQ) 


(4) 


gi^B gi{Si!/>s+Qi<l>Q) 

B 


, g-i^B ^-i{Si4‘S+Qi<pQ) ^ 


where baryonic and antibaryonic terms have been sepa¬ 
rated. Introducing the notation J2b = 


w, the above equation can be rewritten as: 


^^(,Bitt>B+Si<j>s+Qi<l>Q) 

S,B 


— '\(bn * 

■ e 




(5) 


Substituting this expression in Eq. © and changing the 
baryon variable according to 4>b ^ 4>b — argw allows 
us to perform analytically the integration 'm (j)B^ which 
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yields a modified Bessel function I^: 

ry pIlT p‘2'K 

Zb,s,q{T,V) = dcj^s d^Q (6) 

X cos(5(/>s+ Q</)q - 5arga;(0s,</)Q)) 

X /b(2|u;(05,(/»q)|) 


X 


exp 


2 ^ zl cos{Si(l)s + Qi(I)q) 

. M 


Thus, we are left with a double integration only which can 
be then performed numerically with no major problem 
provided that, B being very large, the uniform asymp¬ 
totic expansion of Bessel functions for large orders is used 
to compute 1b- 

It is worth noting, that the form does not involve 
a saddle-point approximation used in 11113113. Ac¬ 
tually, for very small systems it is not even applicable, 
because the integrand © is rather smoothly distributed 
in the group parameter space {4>b,4'Sj4‘q)j so the basic 
assumption of the saddle-point method fails. 

The mean particle numbers of primary hadrons (i.e. 
those directly emitted from the hadronisin g so urce) (Ni) 
are obtained from the equation by taking 

the derivative of the canonical partition function with 
respect to a fictitious fugacity A^: 




, d\nZB,s,Q{T,V) 

ZB-Bi,S-Si,Q-Qi{T, V) 1 
ZB,S,QiT,V) 


(7) 


where the quant^ ZB-Bi,s-Si,Q-QilZB,s,Q is called 
chemical factor [3- In the large volume (thermody¬ 
namical) limit, this expression becomes the GC one, i.e. 
{Ni) = Af ‘ A|* Aq’. 

The canonical baryon chemical factor can be defined 
as that relevant to baryons with vanishing strangeness 
and electric charge, e.g. neutrons: 


Zb-i,s,q{T, V) 
ZB,s,QiT, V) 


( 8 ) 


so that in the GC limit one gets: 


lim Cb = Ab. (9) 

V^OO 


Similarly, one can define chemical factors Cs and Cq for 
strangeness and electric charge respectively. It must be 
emphasized that, unlike in the GC framework where the 
total fugacity is actually a product of different charge fu- 
gacities powered to the number of charges carried by the 
hadron, the canonical chemical factors cannot be factor¬ 
ized and must be calculated, as indicated in equation 0 , 
as a ratio of partition functions for different sets of quan¬ 
tum numbers. The relative differences between chemical 
factors and their GC counterparts and quantities related 
to these differences are defined, in this paper, as canoni¬ 
cal effects. 


In addition to the comparison between full canonical 
calculations and GC ones, it is very important to assess 
whether the enforcement of exact strangeness conserva¬ 
tion with GC treatment of baryon number and electric 
charge {strangeness canonical ensemble) is indeed suf¬ 
ficient to produce results close enough to full canoni¬ 
cal ones. Here, we only quote the analytical form of 
strangeness-canonical partition function Zs derived us¬ 
ing the method introduced in M 

Zs{T,V,Xb,Xq) = (10) 

oo oo 

Zo ^ /„(2 Zho)AS ^ /„(2Zh-)ASAq" 

m— — oo n— — oo 

oo 

X ^ /i(2Zn-)AUQ 
1 — — 00 

S-\-2m-\-2n-\-3l 

Is+2m+2n+3l (2 \/A/lA/’_ i ) . 



In the equation above, we denote the sum of one par¬ 
ticle partition functions carrying the same quantum 
numbers as 'EP by and similarly for other kinds 

of multistrange hadrons, while A/) stands for the sum 
^b'^q' ^i ■ case, mean particle numbers 

read: 


{N^) 


Zs-Si{T,V,XB,XQ) ^B, 1 

Zs{T,V,Xb,Xq) ^ 


( 11 ) 


III. NUMERICAL RESULTS 

First of all, we point out that we do not perform fits to 
experimental results in this paper. Indeed, as our main 
motivation is to test the statistical canonical formalism in 
the heavy ion collision regime, we only take some sugges¬ 
tive values for freeze-out temperatures and baryon den¬ 
sities from our recent work on the analysis of particle 
multiplicities in heavy ion systems in the laboratory mo¬ 
mentum range from GSI SIS 1.7 A GeV to CERN SPS 
158 A GeV [3 and make comparisons between canonical 
and grand-canonical calculations for those values. 

The peculiar feature of canonical formalism is a nonlin¬ 
ear volume behaviour of mean particle numbers because 
of the dependence of chemical factors on the volume ac¬ 
cording to eqs. fTia . To show this effect, we first fix the 
baryon density ub so that the dependence of the chemi¬ 
cal factors on the volume turns into a dependence on net 
baryon number. As the temperature is fixed too, then the 
grand-canonical baryon chemical potential can be calcu¬ 
lated and the difference with the corresponding chemical 
factor can be studied. Furthermore, net strangeness S 
and electric charge Q are fixed by the initial conditions 
of the collision (namely zero and {ZjA) x B respectively). 

The canonical effect is the most significant at compar¬ 
atively low energies, such as 1.7 GeV per nucleon in SIS 
Au-Au collisions. In our recent analysis Q we find the 













4 



FIG. 1: Canonical baryon chemical factor as a function of 
number of participants in the conditions relevant to GSI en¬ 
ergies. Thermodynamical limit is shown as a dashed line. 
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FIG. 2: Canonical enhancement of kaons as a function of num¬ 
ber of participants in the conditions relevant to GSI energies. 


best fit to experimental results with thermal parameters 
T ~ 50 MeV and /rs ~ 800 MeV. These values yield, 
in the pointlike hadron gas framework, a baryon density 
lower than 0.1 fm“^, which is used as a reference value 
in order to study the canonical effect in heavy ion colli¬ 
sions as a function of the number of participants - i.e. net 
baryon number - in the collision. In figure ^ the chemi¬ 
cal factor Cs is shown in comparison with its GC limit 
Xb- Under these circumstances, the canonical effect re¬ 
lated to the baryon chemical factor is found to decrease 
from 2.1% at B = 2 down to 0.7% at B = 10, and to 
further decrease rather slowly towards the GC limit with 
increasing baryon number or volume. For central Au-Au 
reactions, the number of participating nucleons is more 
than 300, so we can conclude that the canonical baryon 
chemical effect is negligible. However, this is not the case 
for strangeness. In figure|21the theoretical production ra¬ 
tio K'*" over number of participating nucleons is shown to 
point out the canonical strangeness effect. Note that we 
have been able to compute the full chemical factor even 
for B = 400. The ratio increases very slowly towards 
the GG limit, thus strangeness must always be handled 
canonically when analysing SIS results. The relative dif¬ 
ference between full canonical and strangeness-canonical 
results is 13% for B = 2, decreasing to 4% for B = 10 
and to 1% for B = 40. It must be pointed out that the 
above results are calculated using an isospin symmetric 
initial configuration whereas Z/A ~ 0.4 for gold nucleus. 
However, this variation essentially gives no change on the 
relative differences above, although the ratio K+/A^p^rt 
decreases naturally due to the initial neutron excess. 

In Au-Au collisions at AGS energies (11.6 A MeV mo¬ 
mentum), the chemical freeze-out temperature is found 
to be around T = 120 MeV and the baryon density close 
to the normal nuclear density Q. Using these values, we 
show the canonical baryon chemical factor in compari¬ 
son with the GG limit in figure 13 keeping a symmetric 
isospin initial condition. Relative difference between Cb 



FIG. 3: Ganonical baryon chemical factor as a function of 
the number of participants in the conditions relevant to AGS 
energies. Thermodynamical limit is shown as a dashed line. 


and Xb decreases rather quickly from 23% at H = 2 to 
1% at i? = 40. A more realistic comparison of results in 
Au-Au with pp collisions requires different initial isospin 
configuration to be into account. By using temperature 
and baryon density quoted above, CB8tB = Q = 2 
turns out to be 81.9 and Xb = 118 with Q/B — 0.4. 

Recent experimental results on AGS K/Aipp,* ratios 
with increasing centrality serve as a basis for a 
further study of strangeness production behaviour at 
T = 120 MeV. In figure^we plot the calculated K+/Vp„t 
with results obtained in peripheral to central Au-Au, Si- 
Au and Si-Al collisions. It can be seen that all the way up 
from H = 2 to H = 60 the full canonical and strangeness- 
canonical results are essentially the same. All central re¬ 
action results lie on the region where the canonical effects 
are negligible and the GC formalism applies. Strangeness 
enhancement in Au-Au system does not definitely look 
like being of canonical origin, whereas Si-Au and Si-Al 
follow roughly the curve at the normal nuclear density 
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FIG. 4: Theoretical K^/A'part curves at fixed temperature for 
two different baryon densities shown along with AGS exper¬ 
imental results IM . Curves are strangeness canonical while 
squares are full canonical results 


FIG. 6: Strange baryon enhancement in the conditions rele¬ 
vant to the SPS energies. Hadron multiplicities are normal¬ 
ized to full canonical results at B = 2. Curves are strangeness 
canonical while crosses are full canonical results. 
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FIG. 5: Theoretical K“/Apart curves at hxed temperature 
for two different baryon densities shown along with AGS ex¬ 
perimental results 12(1 . The curves are strangeness canonical 
while squares are full canonical results 


(not shown in the figure). The above argument also ap¬ 
plies for the K“/Apart enhancement pattern shown in fig¬ 
ure 0 

Whereas the kaon enhancement with increasing cen¬ 
trality seems not to be compatible with a purely canoni¬ 
cal effect, inelastic pp and central p-A and A-A systems 
considered in ref. [g show a different behaviour with in¬ 
creasing volume. Therein, it was found that K/tt ratios 
clearly follow the canonical curves. 

The best fit temperature for the multiplicities mea¬ 
sured by NA49 experiment at SPS in central Pb-Pb re¬ 
actions at 158 A GeV is found to be T ~ 160 MeV 0. 
In figure El we show the strange baryon enhancement by 
using a baryon density 0.3 fm“^. Although the baryon 
density in our analysis 0 was found to be ns ~ 0.2 fm“^ 
this larger value is used in order to probe the largest 


reasonable canonical effect. All results in figure El are 
normalized to the full canonical baryon multiplicities per 
participant at the point Ap^^ = 2. This choice reveals the 
slight difference between the results obtained using the 
S'-canonical approximation and the (B, S, Q)-canonical 
calculation, and forces both canonical methods to reach 
for the same GC limit. An interesting feature here is the 
fact, that the S'-canonical method leads to an overestima¬ 
tion of the canonical effect. Gorresponding curves have 
also been calculated in ref. within the strangeness 
canonical ensemble with some numerical approximations 
not used here. This enhancement picture has been calcu¬ 
lated again with the assumption of an initial symmetric 
isospin configuration. Ghanging the freeze-out conditions 
to T = 160 MeV and ub = 0.17 fm“3, the total relative 
enhancement of (A) from canonical B = Q = 2 to 
GG QjB = 0.4 raises to 30% for A’s, 241% for S’s and 
918% for fl’s. 

SPS NA49 collaboration has also measured the par¬ 
ticipant dependence of K~*~ / 7 r~* ~ ratio in pp, G-G, Si-Si, 
S-S and Pb-Pb reactions . In fig. [ 7 | we show that 
the roughly linear enhancement pattern is far away from 
the purely canonical shape. In the same figure, one can 
see the weak dependence of the K+/ 7 r+ ratio on thermal 
parameters. Indeed, theoretical curves always lie well 
above the measured points, which is a strong indication 
that strangeness is not in complete chemical equilibrium. 
Usually, this lack of equilibrium is taken into account by 
introducing a parameter 75 (see discussion in section llVll , 
which, in the Boltzmann approximation, appears as a 
linear coefficient in front of the K/tt ratios. Taking our 
best fit value 75 ^ 0.8 0 | brings the theoretical curves 
asymptotically to the experimental central Pb-Pb reac¬ 
tion points. Gurves with T = 160 MeV (lowest two) 
assume initial isospin symmetry, whilst the S canonical 
curve with T = 150 MeV is corrected for the condition 
Q/B = 0.4. It is worth quoting our result for pp to cen- 
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FIG. 7: Theoretical curves at fixed temperature 

compared with SPS experimental results Curves are 

strangeness canonical while crosses are full canonical results. 
Curves from top to bottom are indicated in the lower right 
corner. Crosses correspond to the curves with T = 160 MeV. 

tral Pb-Pb K+/ 7 r+ = R enhancement using T = 160 
MeV and ub = 0.17 fm“^ taking into account the initial 
isospin: i?pb_pb/i?pp = 1.47 with 75 = 1. Experimentally, 
this ratio is about 2 as it can be inferred by looking at 
the points in fig. [ 7 | 


IV. DISCUSSION ON THE RELATIVE 
CHEMICAL EQUILIBRIUM 


In statistical model analysis of heavy ion reactions it 
is found that the assumption of full chemical equilibrium 
of strangeness is mostly not satisfied ( 0 ] and this work). 
Thus, a concept of relative chemical equilibrium was in- 
roduced and has been applied in many subsequent 
works. The relative chemical equilibrium means that the 
strange quark production is not fully equilibrated, but 
the distribution of strange quarks among different hadron 
species occurs according to statistical equilibrium {statis¬ 
tical coalescence). 

In order to parametrize the relative chemical equilib¬ 
rium, another constraint in the derivation of the relevant 
statistical operator is needed. In the GC formalism, this 
constraint is the conservation of the average number of 
strange and antistrange quarks in addition to usual con¬ 
straints for energy and charge conservation. Note that 
this constraint is not owing to any natural symmetry of 
strong interactions like baryon number, electric charge 
and strangeness conservation. Hence, as the number of 
produced strange quarks to be distributed among final 
hadrons is unknown, it always leads to the introduction 
of another free parameter to be determined a posteriori. 

Using the aforementioned constraints, the grand 
canonical (GC) partition function can be written as: 


ZGG{T,V,p.,p,n,) = Tr 




( 12 ) 


where the chemical potentials for all the relevant 
charges have been grouped under the symbol p. The 
fugacity factor for the constrained number of strange 
quarks is usually called 75 

In the full canonical picture, where charges are not 
allowed to fluctuate from sample to sample, the con¬ 
straint of the average number of strange and antistrange 
quarks can be turned into the enforcement of a fixed ex¬ 
act number of s, s quarks. Corresponding U (l)™^ symme¬ 
try yields another integration variable and another phase 
factor, exp(ms(/)„^), into the projection integral |(2J) [l7| . 
Following the procedure described in section |n] the aver¬ 
age number of hadrons hi can be written as: 

< hi >- -^. (Id) 

^B,S,Q,ns 

In the grand canonical limit, the corresponding number 
is 


lim < hi >— ^b'^ s' 
V—^OO ^ 


As an illustrative example, take the average number of (j) 
mesons. It does not carry any relevant charge, but the 
valence quark content is ss. In the GC limit <(/)>= 
7520 . This leads us to the identification 


lim 

V^—>^00 


ZB,S,Q,ns-2 

Zb.S.Q.u, 


= Is- 


(15) 


There is one major caveat in those formulae and the 
use of 75 . Indeed, the canonical formalism for the net 
conserved charges in the system is to be applied in any 
reaction where they are known from the initial state; if 
the system is large enough, the fluctuations allowed in 
GC ensemble give negligible deviation from canonical re¬ 
sults and one is allowed to use formulae such as eq. Cl- 
However, the case of the number of strange quarks Ng is 
instrinsically different. In fact, this number is not known 
from the initial state and may undergo large dynami¬ 
cal fluctuations from event to event, well beyond those 
(small) predicted by the GC formalism and its related 
fugacity 75 . If this was the case, 7 s in eq. da must be 
understood as a sort of average fugacity rather than a 
proper fugacity. On the other hand, when the expected 
mean number of ss pairs is small, the canonical formalism 
should be used (i.e. eq. da) but the probability distri¬ 
bution of generating a given number of ss pairs should 
be known in advance. As a reasonable ansatz, one can 
assume those pairs are independently produced so that 
their distribution is a Poissonian [l^ . 


V. CONCLUSIONS 

We have presented a method to calculate particle den¬ 
sities in the canonical framework of the statistical model, 
which can be effectively used in heavy ion reactions at 
very large values of baryon number. This method has 
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allowed us to study the applicability of different approxi¬ 
mations in chemical analysis, namely the grand-canonical 
and the strangeness-canonical ensembles. It is found that 
the full {B, S, Q) canonical formalism is only needed for 
elementary and very small nuclear systems up to baryon 
number 10. For larger systems, the baryon number 
and the electric charge can be safely handled in grand- 
canonical manner. For central Au-Au and Pb-Pb reac¬ 
tions, exact strangeness conservation is needed at GSI 
SIS energies < 2 GeV A), while at AGS and SPS 
even strangeness can be handled grand-canonically. Ger- 
tainly, this also applies to much higher energy RHIG and 
LHG collisions. 

The observed strangeness enhancement from periph¬ 
eral to central nuclear collisions is generally not 

reproduced by the canonical curves with fixed tempera¬ 
ture and baryon density whilst the AGS p-p to p-A to 
A-A central reaction measurements are found to follow 
the canonical results Q. This is an indication of a possi¬ 
ble change in reaction systematics going from peripheral 
to central collisions. The strange hadron enhancement 
in SPS from p-p to Pb-Pb reactions is only qualitatively 


reproduced by canonical model. 

Whilst 7 s was not used in the derivation of our 
numerical results, the weak dependence of K/tt ratio 
on the thermal parameters (shown in fig. [71) serves 
as a further compelling piece of evidence in favour of 
incomplete strangeness chemical eqilibration in GERN 
SPS heavy ion reactions, confirming previous findings 

UM 
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